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Introduction 


A  simulation  model  for  biomass  dynamics  of  a  submersed  Eurasian 
watermilfoil  vegetation  has  been  developed,  called  MILFO.  The  model  is  based 
on  carbon  flow  through  the  vegetation  within  a  1-m^  water  column.  It  includes 
several  aspects  that  affect  biomass  dynamics,  such  as  site-characteristic  changes 
in  climate,  water  temperature,  pH  and  oxygen  effects  on  CO2  assimilation  rate  at 
light  saturation,  wintering  strategies,  grazing  and  mechanical  control  (removal  of 
shoot  biomass),  and  of  latitude.  The  characteristics  of  the  community  and  of  site 
can  be  easily  modified  by  the  user.  MILFO  is  based  on  modeling  concepts  and 
approach  similar  to  those  used  to  model  another  submersed  macrophyte,  hydrilla 
(HYDRIL;  Best  and  Boyd  1996;  Boyd  and  Best  1996). 

MILFO  incorporates  insights  into  the  processes  affecting  dynamics  of  a 
Eurasian  watermilfoil  community  in  relatively  shallow,  hard  water  (0.5-6  m 
depth;  DIG  concentration  >  0.8  mmol  and  pH  ranging  from  7.6  to  9.4).  It  has 
been  calibrated  on  data  pertaining  to  a  milfoil  community  in  Lake  Wingra, 
Wisconsin.  At  that  site,  growth  starts  from  the  basal  rhizome/root  crown  system, 
alone  or  with  wintering  shoot  biomass  present.  Shoot  biomass  usually  peaks 
twice  a  year,  in  June  originating  from  the  first  plant  cohort  and  in  August  from 
the  second  plant  cohort,  and  intensive  downward  transport  of  soluble 
carbohydrates  occurs  after  each  flowering  period.  In  tropical  climates,  a  third 
plant  cohort  is  active. 

MILFO  simulated  the  dynamics  of  plant  and  rhizome/root  crown  biomass  at 
Lake  Wingra  well  over  a  period  of  1  to  5  years.  It  has  been  used  to  calculate  plant 
and  rhizome/root  crown  biomass  for  the  same  latitude  in  a  different  year  and  for 
other  latitudes  in  temperate  (Alabama,  USA)  and  tropical  (India)  areas,  where  it 
gave  good  results. 

Sensitivity  analysis  showed  that  maximum  plant  biomass  of  a  Eurasian 
watermilfoil  community  is  most  sensitive  to  changes  in  light-use  efficiency  and 
very  sensitive  to  photosynthetic  activity  at  light  saturation,  and  that  end-of-year 
rhizome/root  crown  biomass  was  usually  more  sensitive  than  maximum  plant 
biomass.  The  latter  illustrates  the  utmost  importance  of  the  rhizome/root  crown 
system  for  local  survival  and  biomass  production  in  milfoil. 

Environmental  factor  analysis  indicated  that  changes  in  climate  greatly  affect 
simulated  maximum  plant  and  rhizome/root  crown  biomass.  Maximum  plant 
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biomass  proved  sensitive  to  changes  in  water  transparency  and,  to  less  extent,  to 
changes  in  water  depth. 

MILFO  can  be  used  as  a  tool  to  predict  the  dynamics  of  a  Eurasian 
watermilfoil  community  over  1  to  5-year  periods.  Running  the  model  with 
different  parameter  values  specific  for  any  particular  site  and/or  treatment,  for 
example,  biomass  removal  to  a  certain  water  depth,  helps  in  gaining  insight  into 
the  predominant  mechanisms  regulating  submersed  plant  dynamics. 

A  detailed  description  of  the  model  is  given  by  Best  and  Boyd  (in 
preparation). 
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2  Installation  and  Execution 
of  the  Model 


System  Requirements 

There  are  few  requirements  for  running  MILFO  (Version  1 .0).  The 
minimum  RAM  memory  requirement  should  be  at  least  5 1 2  kb.  A  mathematical 
coprocessor  is  in  general  not  required  but  will  often  speed  up  the  calculations 
considerably.  A  free  hard  disk  space  of  about  1  Mb  is  required. 


Installing  MILFO  (Version  1.0) 

Create  a  new  directory,  with  the  name  “MILFOIL,”  using  the  DOS  command 
MKDIR.  For  example,  the  DOS  command  “MKDIR  C:\MILFOIL”  creates  a 
new  directory  on  the  C  drive  called  MILFOIL.  Copy  the  contents  of  the  floppy 
diskette  to  the  directory  C:\MILFOIL  by  using  the  following  command: 

XCOPY  A:\MILFOIL\*.*  C:\MILFOIL\  *.*  /s 

This  diskette  contains  the  MILFO.EXE  file,  all  necessary  input  data  files,  as  well 
as  a  file  used  to  display  the  model  output  graphically  (TTSELECT.EXE).  Input 
files  included  with  this  diskette  are  as  follows: 

a.  MODEL.DAT 

b.  TIMER.DAT 

c.  CONTROL.DAT 

d.  RERUNS.DAT 

Available  weather  data  files  are  included  in  the  subdirectory  AVEATHER. 
The  user  can  select  any  one  of  these  weather  files  as  input  for  the  MILFO  model 
or  may  choose  to  create  a  weather  file  specific  to  a  particular  site.  The  content  of 
the  weather  data  files  is  discussed  in  Chapter  3  of  this  manual,  and  an  example  of 
a  typical  weather  data  file  can  be  seen  in  Appendix  A. 
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Executing  MILFO  (Version  1.0) 

The  MILFO  model  does  not  require  interactive  input  during  execution.  The 
runs  have  been  specified  completely  in  the  data  files.  To  execute  the  model, 
simply  type 


MILFO  <CR>  (carriage  return) 

An  introductory  screen  appears,  and  the  user  is  prompted  to  press  <ENTER>. 
During  execution,  the  model  will  display  the  run  number,  year  number,  and  day 
number  on  the  screen  each  time  output  to  file  is  done.  During  execution,  errors 
and  warnings  may  occur  from  the  weather  system  and/or  from  the  other  modules 
of  the  model.  These  errors/wamings  generally  consist  of  one  line  of  text.  If  the 
simulation  is  terminated  by  an  error  during  the  dynamic  section  of  the  run,  the 
outputs  generated  before  the  error  in  that  particular  run  occurred  are  written  to  a 
temporary  file  but  are  not  written  to  the  output  file  until  the  terminal  section  of 
the  model  is  reached. 
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3  Program  Structure  and 
Data  Files 


Program  Structure 

The  source  code  for  MILFO  (Version  1 .0)  is  written  in  Fortran??.  The 
model  runs  within  a  system  called  the  FORTRAN  Simulation  Environment 
(Version  2.1 ),  hereafter  referred  to  as  FSE.  The  FSE  allows  a  simulation  model 
to  be  written  with  emphasis  on  the  modeling  itself  rather  than  on  other  things 
such  as  time,  file  i/o,  etc.  More  information  is  available  on  running  the  model 
within  the  FSE  in  Chapter  5  of  this  manual. 

Subroutines  called  during  the  execution  of  the  MILFO  model  include 
MODELS,  MODEL,  ASTRO,  TOTASS,  ASSIM,  CHORT2,  and  CHORT3.  A 
brief  description  of  each  subroutine  follows: 

MODELS  -  This  subroutine  is  the  interface  routine  between  the  FSE-driver  and 
the  simulation  model.  The  FSE-driver  calls  this  routine  and  transfers  relevant 
“environment”  variables  (such  as  TIME,  OUTPUT,  etc.)  to  this  routine. 

MODEL  -  This  subroutine  is  called  from  subroutine  MODELS.  It  is  the  routine 
where  specific  calculations  for  milfoil  growth  begin. 

ASTRO  -  This  subroutine  is  called  from  the  MODEL  routine  each  day  of  the 
simulation  period.  It  calculates  astronomic  day  length,  photoperiodic  day  length, 
and  diurnal  radiation  characteristics. 

TOTASS  -  This  subroutine  is  called  in  the  MODEL  subroutine.  This  subroutine 
calculates  daily  total  gross  assimilation  by  performing  a  Gaussian  integration 
over  time.  At  three  different  times  of  the  day,  radiation  is  computed  and  used  to 
determine  assimilation. 

ASSIM  -  This  subroutine  is  called  from  subroutine  TOTASS.  Plant  biomass  is 
distributed  within  the  layers  of  the  plant,  and  the  instantaneous  carbon  dioxide 
assimilation  rate  of  the  plant  is  computed  in  this  subroutine. 

CHORT2  -  This  subroutine  is  called  from  subroutine  MODELS.  It  is  the  routine 
where  specific  calculations  for  milfoil  growth  of  COHORT2  begin. 
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CHORT3  -  This  subroutine  is  called  from  subroutine  MODELS.  It  is  the  routine 
where  specific  calculations  for  milfoil  growth  of  COHORTS  begin. 

A  diagram  illustrating  the  program  structure  of  the  model  is  shown  in  Figure  1. 


Data  Files 

Most  of  the  parameters  and  initial  values  of  the  various  processes  are  read 
from  data  files.  This  has  the  advantage  in  that  the  model  does  not  have  to  be 
recompiled  and  linked  each  time  changes  are  implemented  to  the  input  data. 
There  are  four  input  files  required  to  run  the  MILFO  model  (excluding  the 
weather  data  file)  with  a  potential  for  other  input  files  available.  The  model  also 
typically  creates  three  output  data  files.  The  input  and  output  files  associated 
with  MILFO  are  discussed  in  this  section. 


Input  files 

The  MODEL.DAT  file:  This  data  file  contains  initial  constants,  model 
parameters,  as  well  as  data  used  for  functions.  An  example  of  the  MODEL.DAT 
file  is  shown  in  Appendix  A. 

The  TIMER.DAT  file:  This  data  file  specifies  variables  for  the  following: 
a.  Time  control. 

( 1 )  Start  time  and  finish  time. 

(2)  Time-step  integration. 

(3)  Year. 
h.  Output. 

(1)  Time  between  different  outputs. 

(2)  Format  of  the  output  file. 

(3)  Selection  of  output  variables. 
c.  Weather  control. 

( 1 )  Directory  in  which  the  weather  data  are  stored. 

(2)  Country  code. 

(3)  Station  number. 


6 


Chapter  3  Program  Structure  and  Data  Files 


MAIN 


^Displays  modal  header 


CaUst=SE 


'END 


I 


subroutine 
MODEL 
Dedaration  of: 

-Formal  parameteis 
■(shell  terms) 

•Modal  parameters  etc. 
-CakwIaOons 
-Cals  to  subroutines 


subrouGne 

subrouGne 

subrouGne 

ASTRO 

TOTASS 

ASSIM 

*  Ught  altenuaGon 

Photoperiod 

Daily  total  gross 

^ - 

water  *  vegetation 

assimiiaUon 

*  Daylength 

{3-polnt  Gaussian 

‘instantaneous 

Inlegration) 

gross  asslmGation 

IrradtaUon 

(0.1  mdepGi  layers) 

.  ^ 

Figure  1 .  Relational  diagram  illustrating  the  organization  of  the  model  MILFO  and  its 

subroutines  in  combination  with  the  FSE  shell  (Each  plant  cohort  is  represented  by  a 
cohort-specific  subroutine  (Cohort  1  by  MODEL,  Cohort  2  by  CHRT2,  and  Cohort  3 
by  CHRT3;  only  one  shown),  all  using  the  same  subroutines,  ASTRO,  TOTASS,  and 
ASSIM) 
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An  example  of  the  TIMER.DAT  file  is  shown  in  Appendix  A. 

The  RERUNS.DAT  file:  If  the  reruns  file  is  absent  or  empty,  the  model  will 
execute  a  single  run  (1  year)  from  the  standard  data  files.  By  creating  a  reruns 
file,  the  model  will  execute  additional  runs  with  different  parameters  and/or 
initial  values  for  the  state  variables  (or  even  different  input  files).  Therefore,  the 
total  number  of  runs  made  by  the  model  is  always  one  more  than  the  number  of 
rerun  sets  (see  Appendix  A).  The  format  of  the  rerun  files  is  identical  to  that  of 
the  other  data  files,  except  that  the  names  of  variables  may  appear  in  the  file 
more  than  once. 

The  CONTROL.DAT  file:  This  file  contains  the  names  of  both  input  and 
output  files  used  during  the  execution  of  MILFO.  An  example  of  the 
CONTROL.DAT  file  is  shown  in  Appendix  A. 


Output  files 

MILFO  creates  three  standard  output  files  with  a  potential  fourth,  binary  file: 
RES.DAT,  MODEL.LOG,  WEATHER.LOG,  and  RES.BIN. 

The  RES.DAT  file:  This  file  contains  the  output  of  the  model  with  the  reruns 
(if  present)  merged  below  each  other  in  the  file.  The  output  variables  are  given 
in  Appendix  B.  The  format  of  the  output  file  RES.DAT  depends  on  the  value  of 
the  variable  IPFORM  from  the  timer  file  (see  Appendix  A). 

The  MODEL.LOG  file:  This  file  may  contain  the  messages  from  routines 
used  during  the  simulation.  Messages  about  replacements  by  the  reruns  facility 
can  be  particularly  useful.  To  make  sure  the  execution  of  the  model  is  without 
errors,  this  file  should  be  inspected. 

The  WEATHER.LOG  file:  This  file  contains  all  the  messages  generated  by 
the  weather  system.  By  default,  all  the  comment  headers  of  the  data  files  and  all 
warnings  and  all  errors  from  the  weather  system  are  written  to  this  log  file.  If 
errors  or  warnings  occur  during  a  run,  a  message  is  displayed  shortly  before  the 
termination  of  the  model  about  possible  errors  or  warnings.  These  messages  are 
explained  in  more  detail  in  the  log  file. 

The  RES.BIN  file:  The  variable  “DELTMP”  found  in  the  TIMER.DAT  file 
(Appendix  A)  determines  if  the  temporary  output  data  (RES.BIN)  should  be 
deleted  or  saved  at  termination  of  the  simulation  (DELTMP  =  ‘N’,  Do  not 
Delete,  DELTMP  =  ‘Y’,  Delete).  Using  this  file,  it  is  possible  to  generate  graphs 
of  the  model’s  output  on  IBM-PCs  and  compatibles  after  termination  of  the 
simulation.  This  can  be  done  using  the  TTSELECT  program,  provided 
DELTMP  is  set  to  ‘N’  in  the  timer  file.  This  program  is  included  in  this 
distribution  package.  For  more  details  on  TTSELECT,  see  Displaying  Output. 
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Weather  data  files 


The  weather  data  system  basically  consists  of  two  parts:  the  weather  data 
files  and  a  program  to  retrieve  data  from  those  files  (Van  Kraalingen  et  al.  1991). 
A  single  data  file  can  contain,  at  most,  the  daily  weather  data  of  one  meteorologi¬ 
cal  station  for  1  particular  year.  The  country  name  (abbreviated),  station  number, 
and  year  to  which  the  data  refer  are  reflected  in  the  name  of  the  data  file  (e.g., 
NLl  .983  applies  to  data  from  the  Dutch  (NL)  meteorological  station  in 
Wageningen  (1)  for  the  year  1983). 


Daily  values  are  provided  for  the  following  weather  parameters: 


Name 


Unit 


Global  radiation(daily  total) 
Minimum  air  temperature 
Maximum  air  temperature 
Vapor  pressure 
Wind  speed  (daily  average) 
Rain  (daily  total) 


kJ/m^/d 

degrees  Celsius 
degrees  Celsius 
kPa 
m/s 

mm/day 


The  user  can  create  a  weather  data  file  that  is  unique  to  a  particular  site.  The 
file  consists  of  four  parts:  a  file  header  containing  some  explanatory  text,  one  line 
with  location  parameters  of  the  station,  lines  with  measured  data,  and,  optionally, 
so-called  status  lines  giving  information  on  the  way  missing  data  should  be 
handled  by  the  reading  program  (see  Van  Kraalingen  et  al.  1991).  An  example  of 
a  weather  file  can  be  seen  in  Appendix  A. 


Displaying  Output 

The  program  TTSELECT.EXE  is  included  within  this  distribution  package. 
Execution  of  this  program  allows  the  user  to  graphically  view  output  parameters 
stored  in  the  file  RES.BIN.  To  use  this  feature  of  the  package,  after  termination 
of  the  MILFO  simulation,  type  the  following  command: 

TTSELECT  <CR>  (carriage  return) 

A  list  of  all  possible  output  parameters  will  then  be  displayed  at  the  top  of  the 
computer  screen.  The  user  must  select  two  or  more  of  these  parameters  by 
entering  the  parameter  name  separated  by  a  comma  (NOTE:  Parameter  names 
must  be  entered  exactly  as  they  appear  on  the  screen).  The  first  parameter 
entered  (always  TIME)  will  appear  as  the  x-axis  variable,  while  all  other 
variables  entered  will  be  plotted  along  the  y-axis.  Once  all  output  parameter 
names  are  entered,  the  user  must  follow  instructions  on  the  screen  by  pressing  a 
<CR>  (carriage  return).  The  output  graph  will  then  be  displayed. 

There  are  several  options  available  once  the  graph  is  displayed:  (a)  the  plot 
can  be  saved  as  a  file;  (b)  it  can  be  saved  as  a  screen  dump  file  for  later  printing; 
or  (e)  it  can  be  printed  on  a  Hewlett  Packard  DeskJet  or  LaserJet  printer.  If 
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desired,  another  set  of  parameters  can  be  viewed  by  entering  different  output 
parameters.  At  anytime,  the  user  may  exit  the  TTSELECT.EXE  program  by 
typing  CONTROL  Z  followed  by  a  <CR>  (carriage  return). 
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4  Program  Output 


Example  Runs 

Example  simulations  using  MILFO  (Version  1 .0)  were  made  to  provide 
further  information,  facilitate  proper  execution,  and  to  demonstrate  applications. 
The  following  summary  includes  two  examples:  one  using  MODEL.DAT  values 
shown  in  Appendix  A  (nominal  values  with  total  live  biomass  beginning  at 
100  g  m'^),  and  one  in  which  no  plant  cohort  biomass  is  present  at  the  beginning 
of  the  simulation. 

The  user  can  produce  simulation  results  for  a  specific  scenario  by  modifying 
parameters  in  the  MODEL.DAT  file.  MILFO  can  rapidly  provide  information  on 
the  growth  and  development  rate  of  milfoil  over  a  specified  period. 


Example  1:  Nominal  run 

The  current  MODEL.DAT  file  (see  Appendix  A)  contains  data  required  to 
execute  a  nominal  run. 

In  the  MODEL.DAT  file  under  the  section  “Initial  constants,”  only  the  initial 
live  biomass  of  the  rhizome/root  crown  and  of  the  first  plant  cohort  have  values  > 
0.  This  run  begins  with  a  total  live  biomass  of  100  g  DW  m’^,  equally  divided 
between  rhizome/root  crown  system  (50  g)  and  first  plant  cohort  (50  g)  on  Day  1 
of  the  simulation.  Initial  biomass  corresponds  for  the  rhizome/root  crown  system 
with  field  data  measured  in  Lake  Wingra,  Wisconsin,  in  1977  (Smith  and  Adams 
1986)  and  for  plant  cohort  biomass  with  field  data  measured  in  the  same  lake  in 
1970  (Adams  and  McCracken  1974).  Of  total  plant  live  weight,  excluding  the 
rhizome/root  crown  system,  0.47  is  allocated  to  leaves,  0.47  to  stems,  and  0.06  to 
roots  (Budd,  Lillie,  and  Rasmussen  1995;  Best  and  Boyd  1996).  The  weights  are 
calculated  below: 

IWGRJZ  =  0.50  X  100  =  50.0  !  initial  dry  matter  rhizome/root  crown 

system  (g  DW  m'^) 

IWLGl  =  0.47  X  (0.50  x  100)  =  23.5  !  initial  dry  matter  of  green  leaves 

(g  DW  m'^) 

IWRGl  =  0.06  X  (0.50  x  100)  =  3.0  !  initial  dry  matter  of  live  roots 

(g  DW  m'^) 
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IWSGl  —  0.47  X  (0.50  x  100)  =  23.5  !  initial  dry  matter  of  green  stems 

(g  DW  m'^) 

All  other  initial  biomass  values,  live  or  dead,  are  set  to  zero. 

The  simulation  is  executed  for  1  full  year  using  weather  data  contained  in  the 
file  'WIS 1 .970'.  This  file  contains  weather  data  obtained  from  the  weather 
station  at  the  Madison  Airport  (minimum  and  maximum  air  temperature,  vapor 
pressure,  and  wind  speed)  and  a  weather  station  located  at  Lake  Wingra 
(irradiance),  both  in  Wisconsin  and  in  1970;  the  file  is  located  on  the 
subdirectory  C:\MILFOIL\WEATHER. 

For  this  nominal  run,  the  model  can  be  executed  without  changing  any  input 
files  used  by  the  model.  To  ensure  that  the  current  directory  is  MILFOIL,  from 
the  C  drive,  type  the  command: 

CD/MILFOIL  <CR> 


and  then  type 


MILFO  <CR> 

to  execute  the  simulation  model. 

After  the  introductory  screen  is  displayed  and  execution  begins,  the  model 
will  display  the  run  number,  year  number,  and  day  number  on  the  terminal  screen 
each  time  output  to  file  is  done.  Upon  completion,  MILFO  lists  the  names  and 
contents  of  output  files  created. 

The  calculation  given  above  has  to  be  done  each  time  a  run  with  a  different 
initial  biomass  is  desired.  The  calculated  initial  values  for  plant  organ  weights 
must  be  changed  in  the  MODEL.DAT  file;  the  MODEL.DAT  file  should  be 
saved;  and  MILFO  should  be  executed  as  described  for  the  nominal  run. 


Example  2:  Initial  rhizome/root  crown  biomass  only 

In  Example  2,  a  MODEL.DAT  file  is  used,  almost  identical  to  the 
MODEL.DAT  file  for  the  nominal  run;  only  initial  plant  biomass  values  under 
the  section  “Initial  constants”  have  been  changed.  This  run  begins  with  50  g  DW 
live  biomass  per  m^  on  Day  1  of  the  simulation.  In  this  situation,  growth  starts 
from  the  rhizome/root  crown  system,  and  initial  biomass  of  all  plant  cohorts, 
including  the  first  one,  is  set  to  0. 


Example  3:  Milfoil  populations  at  different  water  depths 

In  Example  3,  a  MODEL.DAT  file  is  used,  almost  identical  to  the 
MODEL.DAT  file  for  the  nominal  run;  only  the  value  for  DEPTH  under  the 
section  “Model  parameters”  has  been  changed  three  times,  i.e.,  for  0.5,  2.5,  and 
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5  m  depth  (nominal  depth  is  1 .5  m).  Changing  this  parameter  enables  running  the 
model  for  the  same  milfoil  stand  in  a  reservoir  with  annually  or  seasonally 
changing  water  depth  and/or  for  different  milfoil  stands  in  lakes  differing  in 
water  depth. 


Example  4:  Milfoil  populations  in  water  bodies  with  different 
transparencies 

In  Example  4,  a  MODEL.DAT  file  is  used,  almost  identical  to  the 
MODEL.DAT  file  for  the  nominal  run;  only  the  values  for  the  function  LT  under 
the  section  “AFGEN  functions”  have  to  be  changed  (results  of  run  not  given). 
Data  pairs  have  to  be  entered  by  giving  first  the  Julian  day  number  followed  by 

and  subsequently  the  value  of  the  extinction  coefficient  at  that  day  followed 
by  Changing  this  parameter  set  enables  running  the  model  for  the  same 
milfoil  stand  in  a  water  for  different  years  with  annually  and/or  seasonally 
changing  light  extinction  coefficient  and/or  for  different  milfoil  stands  in  lakes 
differing  in  a  water  transparency.  Water  transparency  values  expressed  in  Secchi 
depth  (in  m)  should  be  converted  to  light  extinction  coefficients  (in  1/m)  as 
follows:  light  extinction  coefficient  =  1 .65/Secchi  depth.  This  conversion  factor 
has  been  found  to  be  reliable  in  a  Secchi  depth  range  from  0.5  to  2  m  (U.S.  Envi^ 
ronmental  Protection  Agency  1992). 


Example  Output 

The  output  of  the  example  Runs  1  and  2  is  presented  in  Figure  2.  Total  plant 
biomass  of  milfoil  was  usually  higher  for  the  nominal  run  than  for  the  run 
without  initial  plant  cohort  biomass  present.  Live  plant  biomass  peaked  at 
Day  240  with  266  g  DW  m’^  in  the  nominal  run  and  with  149  g  DW  m'^  in  the 
run  without  Initial  plant  biomass  present.  End-of-year  dead  plant  biomass  was 
also  considerably  higher  in  the  nominal  run  (590  g  DW  m'^)  versus  the  run 
without  initial  plant  biomass  present  (317  g  DW  m'^).  Initial  biomass  influenced 
not  only  the  total  live  and  total  (live+dead)  diy  weight  to  a  great  extent  but  also 
influenced  rhizome/root  erown  weight. 

All  output  parameters  listed  in  Appendix  B  can  be  displayed  by  selecting  the 
desired  abbreviation  using  TTSELECT,  since  all  data  are  contained  in  the 
RES.BIN  fde.  However,  if  the  numbers  on,  for  example,  the  accumulated  live 
biomass  of  the  three  plant  cohorts  are  desired  as  output,  the  abbreviation  of  this 
term  (TGW)  should  be  ‘freed’  (remove  star  in  front  of  this  term)  in  the 
TIMER.DAT  file,  and  the  user  can  find  the  desired  data  in  the  RES.DAT  file.  For 
the  accumulated  live+dead  biomass,  the  term  would  be  TW;  for  the  live  weights 
of  the  plant  organs  of  one  individual  plant  cohort,  the  terms  would  be  TWLGl, 
TWSTl,  TWRGl,  and  for  their  dead  weights,  TWLDl,  TWSDl,  TWRDl. 

The  output  of  example  Run  3  is  presented  in  Figure  3.  Simulated  live  plant 
and  rhizome/root  crown  biomass  generally  decreased  with  increasing  water 
(=rooting)  depth. 
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We  hope  that  these  example  runs  and  their  output  provided  in  this  section  of  the 
users’  manual  will  help  users  conceptualize  simulation  runs  that  will  generate 
useful  information  to  address  specific  questions. 


Use  of  Different  Weather  Files 

Climate  is  an  important  factor  influencing  milfoil  biomass.  Phenology  is  tied 
indirectly  to  air  and/or  water  temperature  through  development  rate.  Weather 
data  are  available  for  different  climatological  conditions  ranging  from  temperate 
to  tropical.  To  illustrate  the  impact  of  climate  on  biomass  of  rhizome/root  crown 
system  and  plants  of  milfoil,  MILFO  was  executed  using  a  weather  file  taken 
from  a  temperate  as  well  as  a  tropical  climate.  In  both  cases,  the  nominal 
MODEL.DAT  file  was  used  for  biomass-specific  input  data.  The  TIMER.DAT 
file  specified  for  the  nominal  run  is  the  temperate  weather  file  (WISl  970)  and 
for  the  tropical  run  the  tropical  file  (INDl  978).  Each  weather  file  is  described 
below: 

Temperate  Weather  File 

Country:  United  States  of  America 

Station:  Madison 

Year:  1970 

Longitude:  89°  33'  W 

Latitude:  43°  13' N 

Elevation:  7  m 

Tropical  Weather  File 
Country:  India 
Station:  Patancheru 
Year:  1978 
Longitude:  78°  28'  E 
Latitude:  17°  27' N 
Elevation:  21  m 

Output  of  both  simulations  is  shown  in  Figure  4,  where  simulated  live  plant 
and  rhizome/root  crown  biomass  in  a  temperate  and  tropical  climate  over  a  1  - 
year  period  is  illustrated.  Apparently,  milfoil  communities  in  a  tropical  climate 
exhibit  a  far  higher  peak  plant  biomass  than  in  a  temperate  climate,  and 
rhizome/root  crown  biomass  tends  to  accumulate  towards  the  end  of  the  growing 
season  in  tropical  regions,  while  it  fluctuates  relatively  more  and  does  not 
accumulate  in  temperate  regions. 


Removing  Plant  Biomass  By  Mechanical 
Harvesting 

MILFO  can  also  be  used  to  calculate  effects  of  various  control  methods  on 
biomass  and  survival  of  a  milfoil  stand,  for  example,  of  mechanical  harvesting  at 
various  times  and  water  depths.  The  model  can  be  run  for  this  purpose  by 
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Temperate:  Lake  Wingra,WI  (1970)  Temperate:  Lake  Wingra.WI  (197 

Tropical:  Patancheru,  India  (1978)  Tropical:  Patancheru,  India  (1978) 
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Figure  4.  Simulated  live  biomass  of  plants  and  rhizomes/root  crowns  of  a  milfoil  community  growing  in  a  temperate  or  a 

tropical  climate  (Initial  biomass  vaiues:  temperate,  shoots,  cf.  Adams  and  McCracken  1974;  rhizomes/root  crowns, 
cf.  Smith  and  Adams  1986;  tropical,  cf.  Zutschi  and  Vass  1973.  Temperate  climate;  data  Lake  Wingra, 
Wisconsin,  1970  (longitude  89°  20’  W,  latitude  43°  08’  N).  Tropical  climate;  data  Patancheru,  India,  1978 
(longitude  78°  28’  E,  latitude  17°  26’  N)) 


making  changes  in  the  section  ‘Model  parameters’  of  the  MODEL.DAT  file,  by 
indicating  (a)  that  harvesting  occurs  (HAR  =  1),  (b)  the  day  at  which  harvesting 
occurs  (HARD AY  =  desired  day  number,  e.g.,  31  July  or  Day  No.  212),  and 
(c)  the  harvesting  depth  below  the  water  surface  (HARDEP  =  desired  depth,  e.g., 
0.8  m).  By  these  changes,  all  plant  biomass  contained  in  those  layers  affected  by 
the  harvesting  depth  is  removed  at  the  specified  day  in  the  simulation  run. 
Examples  of  effects  of  various  harvesting  times  and  depths  are  presented  by  Best 
and  Boyd  (in  preparation). 


Multiple-Year  Runs 

MILFO  has  the  ability  to  generate  multiple-year  simulation  runs.  This  is  a 
critical  feature  when  examining  plant  growth  for  consecutive  years  and/or  in  the 
analysis  of  the  effect  of  a  different  value  for  an  input  parameter.  Multiple-year 
simulation  runs  can  be  accomplished  using  the  RERUNS.DAT  file,  which  is 
illustrated  in  Appendix  A  of  this  manual.  If  the  RERUNS.DAT  file  is  absent  or 
empty,  the  model  will  execute  a  single  run,  using  the  data  from  the  standard  data 
files  (i.e.,  MODEL.DAT  and  TIMER.DAT).  If  the  RERUNS.DAT  file  is  present 
and  contains  different  parameters  and/or  initial  values  for  the  state  variables,  the 
total  number  of  runs  made  by  the  model  is  always  one  more  than  the  number  of 
rerun  sets.  Names  of  variables  originating  from  different  data  files  can  be 
redefined  in  the  same  rerun  file.  Arrays  can  also  be  redefined  in  a  rerun  file.  The 
order  and  number  of  variables  should  be  the  same  in  each  set.  A  new  set  starts 
when  the  first  variable  is  repeated  (see  Appendix  A). 
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5  Running  the  Model  Within 
a  Shell 


This  chapter  gives  a  brief  description  of  how  the  FSE-shell  drives  the 
MILFO  model.  All  execution  starts  with  a  MAIN  program  (Figure  1).  This  is  a 
short  program  that  displays  the  header  and  calls  the  FSE-driver.  The  FSE-driver 
then  performs  a  number  of  actions.  It  reads  the  input  and  output  file  names 
needed  by  the  model  from  the  file  CONTROL.DAT.  This  file  contains  the 
names  of  the  input  files  TIMER.DAT,  RERUNS.DAT,  and  MODEL.DAT.  The 
CONTROL.DAT  file  also  contains  names  of  the  model’s  output  files  (RES.DAT 
and  MODEL.LOG).  From  the  weather  control  variables  in  the  TIMER.DAT  file, 
the  weather  system  determines  which  weather  data  file  is  required. 

The  FSE-driver  then  calls  a  MODELS  subroutine  and  transfers  all  relevant 
“environment”  variables  (such  as  TIME,  OUTPUT,  etc.)  to  this  routine.  The 
MODELS  subroutine  provides  the  interface  between  the  FSE-driver  and  the 
simulation  model.  This  routine  in  turn  calls  the  MODEL  subroutine  that  begins 
execution  of  the  various  routines  within  the  MILFO  source  code. 

It  is  not  necessary  to  know  the  FORTRAN  details  of  what  is  going  on  in  the 
FSE-driver.  A  discussion  of  information  that  is  passed  from  the  FSE-driver  to  the 
model  and  vice  versa  can  be  found  in  Van  Kraalingen  (1995). 
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6  Model  Features 


Features  of  MILFO  are  as  follows: 

a.  Phenology  is  tied  indirectly  to  air  temperature  through  development  rate 
and  is,  therefore,  independent  of  day  of  year;  thus,  the  model  can  be  used 
under  climatological  conditions  ranging  from  temperate  to  tropical. 

b.  Plant  growth  starts  from  the  rhizome/root  crown  system  alone  or  from  the 
same  system  with  wintering  plants. 

c.  Two  plant  cohorts  are  active  in  a  temperate  climate  and  three  cohorts  in 
the  tropics. 

d.  Photosynthetic  response  is  to  instantaneous  irradiance. 

e.  Removal  of  biomass  through  harvesting  can  be  calculated  if  desired. 

/  Air  and/or  water  temperatures  can  be  used  to  run  the  model. 

g.  The  model  can  be  used  for  communities  at  various  water  depths,  ranging 
from  0.5  to  6.0  m. 

h.  Plant  parameter  values  and  climatological  variables  can  be  easily 
changed. 
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7  Application  Possibilities 


MILFO  can  be  used  to  assess  behavior  of  a  milfoil  vegetation  under  various 
site-specific  and  climatological  conditions,  and  it  can  be  run  with  user-specified 
input  values  for  plant  and  rhizome/root  crown  biomass. 

Effects  of  man-made  control  activities  like  harvesting  at  different  times  and 
at  various  water  depths  can  be  calculated  also.  Thus,  in  the  latter  case,  it  can  be 
used  as  a  tool  for  aquatic  plant  management  agencies. 

The  present  version  of  MILFO  (1.0)  has  been  developed  as  a  stand-alone 
simulation  model.  It  can  be  relatively  easily  modified  to  communicate  with 
ecosystem  models  because  it  is  written  in  FORTRAN??  and  its  structure  is 
simple.  It  is  planned  to  link  MILFO  to  a  Geographical  Information  System 
through  an  appropriate  interface  like  AEGIS+  (Luyten  et  al.  1994). 
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Appendix  A 
Examples  of  Data  Files 


Example  of  MODEL.DAT  File 

*  MODEL.DAT  file 

*  contains:  -  Initial  constants  as  far  as  specified  with  INCON  statements, 

*  -  Model  parameters, 

*  -  AFGEN  functions, 

*  -  A  SCALE  array  in  case  of  a  general  translation 

* 

*  File:  MDELMW70.DAT;  to  be  used  as  input  file  for  MILFO.FOR 

*  Contains  data  on  biomass  and  water  temperatures  from  Lake  Wingra,  Wl,  1970 

*  Date:  10-10-97 

*  Time:  10:00:06 

* 

*  Initial  constants 

*  _ 

CRRIZ  =12.  !  Critical  weight  of  the  rhizome/root  crown  system  (g  DW  m'^) 

IWGRIZ  =  50.  !  Initial  weight  of  live  rhizome/root  crown  system  (g  DW  m‘^) 

IWRIZD  =  0.  !  Initial  weight  of  dead  rhizome/root  crown  system  (g  DW  m"^) 

IWLDl  =  0.  !  Initial  dry  matter  of  dead  leaves  plant  cohort  1  (g  DW  m'^) 

IWLD2  =  0.  !  Initial  dry  matter  of  dead  leaves  plant  cohort  2  (g  DW  m'^) 

IWLD3  =  0.  !  Initial  dry  matter  of  dead  leaves  plant  cohort  3  (g  DW  m'^) 

IWLGl  =23.5  !  Initial  dry  matter  of  green  leaves  plant  cohort  1  (g  DW  m'^) 
IWLG2  =  0.  !  Initial  dry  matter  of  green  leaves  plant  cohort  2  (g  DW  m'^) 

IWLG3  =  0.  !  Initial  dry  matter  of  green  leaves  plant  cohort  3  (g  DW  m'^) 

IWRDl  =  0.  !  Initial  dry  matter  of  dead  roots  plant  cohort  1  (g  DW  m"^) 

IWRD2  =  0.  !  Initial  dry  matter  of  dead  roots  plant  cohort  2  (g  DW  m'^) 

IWRD3  =  0.  !  Initial  dry  matter  of  dead  roots  plant  cohort  3  (g  DW  m'^) 

IWRGl  =3.  !  Initial  dry  matter  of  live  roots  plant  cohort  1  (g  DW  m"^) 

IWRG2  =0.  !  Initial  dry  matter  of  live  roots  plant  cohort  2  (  g  DW  m'^) 

IWRG3  =0.  !  Initial  dry  matter  of  live  roots  plant  cohort  3  (  gDW  m'^) 

IWSDl  =  0.  !  Initial  dry  matter  of  dead  stems  plant  cohort  1  (  g  DW  m'^) 

IWSD2  =  0.  !  Initial  dry  matter  of  dead  stems  plant  cohort  2  (g  DW  m'^) 
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IWSD3 

=  0. 

!  Initial  dry  matter  of  dead  stems  plant  cohort  3  (g  DW  m'^ 

IWSGl 

=  23.5 

!  Initial  diy  matter  of  live  stems  plant  cohort  1  (g  DW  m'^) 

IWSG2 

=  0. 

!  Initial  dry  matter  of  live  stems  plant  cohort  2  (g  DW  m"^) 

IWSG3 

=  0. 

!  Initial  dry  matter  of  live  stems  plant  cohort  3  (g  DW  m'^) 

NUL 

=  0. 

!  Zero  (0) 

REMOB  1 

=  0. 

!  Remobilization  rate  of  carbohydrates  plant  cohort  1 

(g  CH2O  m'^  d') 

*  Model  parameters 

*  _ 


YRNUM 

=  1. 

!  Year  number  simulation  (1  -  5)(y) 

AMX 

=  0.0165 

!  Potential  CO2  assimilation  rate  at  light  saturation 
!  for  shoot  tips 
!  (g  CO2  g  DW'  h'') 

CVT 

=  1.05 

!  Conversion  factor  for  translocated  dry  matter 
!  into  CH2O  (-) 

DAYEM 

=  1. 

!  First  Julian  day  number  (d) 

DELAY 

=  7. 

!  Lag  period  chosen  to  relate  water  temperature  to  air 
!  temp.,  in 

!  cases  where  water  temp,  has  not  been  measured  (d) 

DEPTH 

=  1.5 

!  Water  depth  (m) 

EE 

=  0.000011 

!  Initial  light  use  efficiency  for  shoots  (g  CO2  J  ') 

HAR 

=  0. 

!  Harvesting  (0=no  harvesting;  l=harvesting) 

HARDAY 

=  212. 

!  Harvesting  day  number  (d) 

HARDEP 

=  0.8 

!  Harvesting  depth  (measured  from  water  surface;  m) 

NPL 

=  11. 

!  Plant  density  (plants  m'^) 

RC 

=  0.06 

!  Reflection  coefficient  of  irradiation  at  water  surface 
!  (relative;  -) 

RCSHST 

=  12.0 

!  Relation  coefficient  rhizome/root  crown-stem  length 
!  (m  g  DW') 

RDRIZ 

=  0.00042 

!  Relative  death  rate  of  rhizome/root  crown  system 
!  (on  dry  weight  basis;  d’’) 

REDAM 

=  0.5 

!  Reduction  factor  to  relate  AMX  to  pH  and  oxygen 
!  levels  of  the  water  (relative;  -) 

ROC 

=  0.0576 

!  Relative  conversion  of  rhizome/root  crown  into  plant 
!  material  (g  CH2O  g  DW' d  ') 

TBASE 

=  3. 

!  Base  temperature  for  juvenile  plant  growth  (°C) 

TL 

=  0.1 

!  Thickness  each  plant  layer  (m) 

TRAFAC  = 

=  0.35 

!  Translocation  factor  (relative;  -) 

*  AFGEN  functions 

* 

!  Developmental  stage  effect  on  AMX  (relative;  -;  presently  not  active) 

*AMDVST  =  0.001,  1. 

3.5,  1.,  5. 0,1. 

!  Daytime  temperature  effect  on  AMX  (relative;  -) 

AMTMPT  = -30.,  0.00001,  0.,  0.00001,  5.,  0.18,  10.,  0.23, 

15.,  0.40,  20.,  0.63,  25.,  0.78,  30.,  0.95, 

35.,  1.0,  40.,  0.78,  45.,  0.38,  50.,  0.00001 
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!  Dry  matter  allocation  to  each  plant  layer  as  function  of  water  depth  (relative;  -) 
DMPCT  =  1.0,  .10,  2.0,  .16,  3.0,  .17,  4.0,  .10,  5.0,  .08 

!  Development  rate  pre-anthesis  as  function  of  daily  average  temperature;  active 
in  source  code; 

!  written  in  MODEL.DAT  file  as  a  reminder) 

*DVRVT  =  -15.,  0.,  0.,  0.,  30.,  0.022 

!  Development  rate  post-anthesis  as  function  of  daily  average  temperature;  active 
in  source  code; 

!  written  in  MODEL.DAT  file  as  a  reminder) 

*DVRRT  =  -15.,  0.,  0.,0.,  30.,  0.015 

!  Leaf  dry  matter  allocation  to  each  layer  of  the  plant  as  function  of  DVS 
(relative;  -) 

FLT  =  0.,  0.50,  2.3,  0.50,  6.0,  0.50 

!  Fraction  of  total  dry  matter  increase  allocated  to  leaves  as  function  of  DVS 
(relative;  -) 

FLVT  =  0.,  0.47,  2.3,  0.47,  6.0,  0.47 

!  Fraction  of  total  dry  matter  increase  allocated  to  stems  as  function  of  DVS 
(relative;  -) 

FSTT  =  0.,  0.47,  2.3,  0.47,  6.0,  0.47 

!  Fraction  of  total  dry  matter  increase  allocated  to  roots  as  function  of  DVS 
(relative;  -) 

FRTT  =  0.,  0.06,  2.3,  0.06,  6.0,  0.06 

!  Plant  species  specific  light  extinction  coefficient  as  function  of  DVS 
(m'gDW',-) 

KT  =  0.,  0.006,  3.5,  0.006,  6.0,  0.006 

!  Water  type  specific  light  extinction  coefficient  (m  ’) 

LT  =  0.,  1.15,  101.,  1.6,  117.,  1.4,  131.,  1.85,  156.,  1.85, 

173.,  1.8,  187.,  1.8,  215.,  1.8,  243.,  2.0,  257.,  1.4, 

271.,  1.6,  299.,  1.4,  327.,  1.15,  365.,  1.15 

!  Relative  death  rate  of  leaves  as  function  of  DAVTMP  (d"') 

RDRT  =  -30.,  0.042,  0.,  0.042,  50.,0.042 

!  Relative  death  rate  of  stems  and  roots  as  function  of  DAVTMP  (d"') 

RDST  =  -30.,  0.042,  0.,  0.042,  50.,  0.042 

!  Reduction  factor  for  AMX  to  account  for  senescence  of  plant  parts  over 
vertical  axis  of  vegetation  as  function  of  DVS  (relative;  -) 

REDFT  =  0.0,  1.0,  1.0,  1.0,  6.0,  1.0 
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!  Factor  accounting  for  effect  of  daily  average  daytime  temperature  on 

maintenance  respiration 

!(-/C) 

TEFFT  =  -30.,  0.0001,  0.,  0.0001,  5.,  1.0,  20.,1.0,  30.,  2.0, 

40.,  4.0,  45.,  8.0,  50.,  12.0 

!  Total  live  plant  diy  weight  measured  as  function  of  day  number  (g  DW  m  ‘) 
TGWMT  =  1.,  50.,  141.,  50.,  162.,  120.,  172.,  106.,  192.,  105., 

202.,  130.,  223.,  160.,  233.,  180.,  254.,  220.,  264.,  150., 

365.,  50.7 


!  Daily  water  temperature  as  function  of  day  number  (°C,  d) 
WTMPT  =  1.,  3.5,  60.,  3.5,  100.,  9.0,  150.,  22.0,  190.,  25.0, 

220.,  25.0,  250.,  19.0,  300.,  9.0,  340.,  1.6,  365.,  1.6 


Example  of  TIMER.DAT  file 


*  TIMER  file  contains: 

* 

*  -  The  used  DRIVER  and  TRACE  in  case  of  GENERAL  translation 

*  -  The  TIMER  variables  used  in  both  translation  modes 

*  -  Additional  TIMER  variables  in  case  of  GENERAL  translation 

*  -  The  WEATHER  control  variables  if  weather  data  are  used 

*  -  Miscellaneous  FSE  variables  in  case  of  FSE  translation 

* 

*  File:  MILFO.FOR 

*  Date:  10-10-97 

*  Time:  10:55:06 


*  TIMER  variables  used  in  GENERAL  and  FSE  translation  modes 

*  _ 


STTIME=  1. 
FINTIM  =  365 
DELT  =1. 
PRDEL  =  1. 
IPFORM  =  4 


!  start  time 
!  finish  time 

!  time  step  (for  Runge-Kutta  first  guess) 

!  output  time  step 
!  code  for  output  table  format: 

!  4  =  spaces  between  columns 
!  5  =  TAB'S  between  columns  (spreadsheet  output) 
!  6  =  two  column  output 


!  The  string  array  PRSEL  contains  the  output  variables  for  which 
!  formatted  tables  have  to  be  made.  One  or  more  times  there  is  a 
!  series  of  variable  names  terminated  by  the  word  <TABLE>. 

!  The  translator  writes  the  variables  in  each  PRINT  statement  to 
!  PRSEL  =  !  a  separate  table. 
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'DAVTMP', 

'DAYL 

'DDTMP 

'DTEFF 

'DTGA 

'DVS 

TGROS 

'GPHOT', 

'IRS 

'MAINRT', 
'MAIN! 
'REMOB  r, 
'REMOB2', 
'REMOB3', 
'TGW, 
'TGWM', 
'TMPSUM', 
'TRANS  1', 
'TRANS2', 
'TRANS3', 
'TW  ', 
'TWGRIZ', 
'TWLDl ', 
'TWLD2 ', 
'TWLD3 ', 
'TWLGl  ', 
'TWLG2 
'TWLG3  ', 
'TWLVD ', 
'TWLVG ', 
'TWRDl  ', 
’TWRD2 ', 
'TWRD3  ', 
'TWRGl  ', 
'TWRG2 ', 
'TWRG3  ', 
'TWRIZ ', 
'TWRIZD', 
'TWRTD ', 
'TWRTG ', 
'TWSDl  ', 
'TWSD2 ', 
'TWSD3 
'TWSGl  ', 
'TWSG2 ', 
'TWSG3  ', 
'TWSTD ', 
'TWSTG ', 
'WTMP ', 
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'<TABLE>’ 

COPINF  =  "N'  !  Switch  variable  whether  to  eopy  the  input  files 

!  to  the  output  file  CN'  =  do  not  copy, 

!  'Y'  =  copy) 

DELTMP  =  "N'  !  Switch  variable  that  should  be  done  with  the 

!  temporary  output  file  CN'  =  do  not  delete, 

!  'Y'  =  delete) 

IFLAG  =1101  !  Indicates  where  weather  error  and  warnings 

!  go  (1 101  means  errors  and  warnings  to  log 
!  file,  errors  to  screen,  see  FSE  manual) 

*IOBSD  =  1 99 1 , 1 82  !  List  of  observation  data  for  which  output  is 

!  required.  The  list  should  consist  of  pairs 
!  <year>,<day>  combination 

*  WEATHER  control  variables 


WTRDIR  =  'C;\MILFOIL\WEATHERV 
CNTR  =  'WIS'  !  Country  code 
ISTN  =  1  !  Station  code 

lYEAR  =  1970  !  Year 


Example  of  RERUNS.DAT  file 

*  RERUNS  file  ...  to  produce  multiple  runs  ... 

* 

*  File:  MILFO.FOR 

*  Date:  09-12-95 
*Time:  11:00:00 


*  RERUNS  variables  used  in  GENERAL  and  FSE  translation  modes 


* 


YRNUM 

=  2. 

!  2"**  year  of  simulation 

HAR 

=  1. 

!  1  =YES,  harvesting  does  occur 

HARDAY 

=  74. 

!  Harvesting  occurs  on  Julian  day  number  74 

HARDEP 

=  1.0 

!  Harvesting  depth  is  1.0  m 

YRNUM 

=  3. 

!  3'^'’  year  of  simulation 

FLAR 

=  1. 

!  1=YES,  harvesting  does  occur 

HARDAY 

=  181. 

!  Harvesting  occurs  on  Julian  day  number  1 8 1 

HARDEP 

=  0.5 

!  Harvesting  depth  is  0.5  m 
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*  WEATHER  control  variables 


WTRDIR  =  'C:\MILFOIL\WEATHERV 
CNTR  =  'WIS'  !  Country  code 

ISTN  =  1  !  Station  code 

lYEAR  =  1972  !  Year 


Example  of  CONTROL.DAT  file 

*  CONTROL.DAT  data  file  contains: 

*  -  File  names  to  be  used  by  FSE  2.1  * 

*  -  The  input  files  (except  FILEIR)  may  be  used  in  reruns;  up  to 

*  five  input  data  files  may  be  used  (FILEI 1  -5) 


FILEON  ='RES.DAT' 

!  Normal  output  file 

FILEOL  ='MODEL.LOG' 

!  Log  file 

FILEIR  =  'RERUNS.DAT' 

!  Reruns  file 

FILEIT  =  'TIMER.DAT’ 

!  File  with  timer  data 

FILEI  1  ='MODEL.DAr 

!  First  input  data  file 

*  FILEI2  =  ’  ’ 

!  Second  input  data  file  (not  used) 

*  FILED  =  ’  ’ 

!  Third  input  data  file  (not  used) 

*  FILEI4  =  ’  ’ 

!  Fourth  input  data  file  (not  used) 

*  FILED  = ' ' 

!  Fifth  input  data  file  (not  used) 

Example  of  a  weather  data  file 

_ _ * 

*  Country;  United  States  of  America 

*  Station:  Madison^WSO  Airport 

*  Year;  1970 

*  Source:  Wisconsin  State  Climatology 

*  Author:  Pam  Knox 

*  Longitude;  89  20  W 

*  Latitude:  43  08  N 

*  Elevation:  858  m. 

*  WMO-code:-  Phone  608-263-2374 

*  Comments:  Data  obtained  from  Pam  Knox  ..  Converted  by  Will  Boyd  (1996) 

*  Irradiation  trom  John  Titus  (pers.comm.1996;  same  data  as  used 

*  for  WEED). 

*  Irradiation  in  Langley/d  multiplied  *  41 .88  to  get  J/m2/day 

*  Columns: 

*  ======= 

*  station  number 

*  year 
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*  day 

*  irradiation  (kJ  m'^  d  ’) 

*  minimum  temperature  (degrees  Celsius) 

*  maximum  temperature  (degrees  Celsius) 

*  vapor  pressure  (kPa) 

*  mean  wind  speed  (m  s"') 

*  precipitation  (mm  d"') 

*  _ if 


89.33 

43.13 

858. 

0.00 

0.00 

1 

1970 

1 

7413. 

-13.3 

-3.3 

3.007 

1.4 

.0 

1 

1970 

2 

8125. 

-15.0 

-6.7 

2.350 

1.1 

.0 

1 

1970 

3 

9088. 

-18.3 

-8.9 

1.827 

3.9 

.0 

1 

1970 

4 

8376. 

-21.7 

-5.6 

1.649 

2.9 

.0 

1 

1970 

5 

10470. 

-24.4 

-15.0 

0.828 

2.7 

.0 

1 

1970 

6 

9674. 

-26.1 

-16.7 

0.562 

2.4 

.0 

1 

1970 

7 

8208. 

-23.9 

16.1 

0.703 

5.1 

.0 

1 

1970 

8 

9381. 

-24.4 

-18.9 

0.629 

7.3 

.0 

1 

1970 

9 

9758. 

-20.0 

-10.6 

1.028 

4.3 

.0 

1 

1970 

10 

8418. 

-25.0 

-8.9 

0.975 

2.4 

.0 

1 

1970 

11 

4314. 

1 

00 

-5.6 

2.727 

3.0 

.0 

1 

1970 

12 

10512. 

-17.8 

-4.4 

1.922 

3.7 

.0 

1 

1970 

13 

7203. 

-21.7 

-10.0 

1.206 

1.5 

2.5 

1 

1970 

14 

6282. 

-14.4 

-5.0 

2.022 

4.3 

.0 

1 

1970 

15 

6743. 

-5.6 

-2.2 

3.478 

4.1 

.0 

1 

1970 

363 

7455. 

-22.2 

-5.0 

1.827 

2.6 

.0 

1 

1970 

364 

6575. 

-14.4 

-3.9 

2.350 

4.4 

.0 

1 

1970 

365 

4439. 

-8.3 

-1.1 

3.478 

3.3 

.0 

*  Note:  Longitude  and  latitude  in  the  header  are  listed  in  degrees  and  minutes 
while  they  are  listed  in  degrees  only  on  the  first  line  of  the  weather  data. 
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Output  Parameters  Available 


Abbreviation 


Explanation 


Dimension 


DAVTMP  Daily  average  temperature 

DAYL  Day  length 

DDTMP  Daily  average  daytime  temperature 

DTEFF  Daily  effective  temperature 

DTGA  Daily  total  gross  CO2  assimilation  of  the  plant 

DVS  Development  phase  of  the  plant 

EE  Initial  light-use  efficiency  for  shoots 

FGROS  Instantaneous  CO2  assimilation  rate  of  the  plant 

GPHOT  Daily  total  gross  assimilation  rate  of  the  community 

IRS  Total  irradiance  just  under  the  water  surface 

MAINT  Maintenance  respiration  rate  of  the  plant 

MAINRT  Maintenance  respiration  rate  of  the  rhizome/root 

crown  system  vegetation  (relative) 

REMOB  1 ,2,3  Remobilization  rate  of  carbohydrates  cohort  1 ,2,3 
temperature  on  maintenance  respiration 
TGW  Total  live  plant  diy  weight  (excluding  rhizome/root 

crown  system) 

TGWM  Total  live  plant  dry  weight  measured  (field  site) 

TMPSUM  Temperature  sum  after  1  January 

TRANS  1,2,3  Translocation  rate  of  carbohydrates  cohort  1,2,3 
TREMOB  Total  remobilization 

TW  Total  live  +  dead  plant  diy  weight  (excluding  rhizome/ 

root  crown  system) 

TWGRIZ  Total  live  rhizome/root  crown  dry  weight 

TWLD1,2,3  Total  dead  leaf  dry  weight  cohort  1,2,3 
TWLG 1 ,2,3  Total  live  leaf  dry  weight  cohort  1 ,2,3 

TWLVD  Total  dry  weight  of  dead  leaves  2  or  3  cohorts 
TWLVG  Total  dry  weight  of  live  leaves  2  or  3  cohorts 
TWRDl ,2,3  Total  dead  root  dry  weight  cohort  1 ,2,3 
TWRG1,2,3  Total  live  root  diy  weight  cohort  1,2,3 
TWRIZ  Total  live  +  dead  rhizome/root  crown  system  weight 

TWRIZD  Total  dead  rhizome/root  crown  system  weight 
TWRTD  Total  dry  weight  of  dead  roots  2  or  3  cohorts 


°C 

h 

“C 

°C 

g  C02.m-ld'’ 

gC02.J-' 
g  C02.m'^h'’ 
g  GHjO-m-ld’’ 
J.m'^.s"’ 
g  CH20.m'^-d-’ 
g  CHjO.m'^-d'’ 

g  CH20.mld’ 

g  DW-m’^ 

g  DW.m’^ 

°C 

g  CH20.mld'’ 
g  CH20.m'^ 
g  DW.m'^ 

gDW.m'^ 
g  DW.m' 
g  DW-m'^ 
g  DW.m'^ 
g  DW.m-^ 
g  DW.m'^ 
g  DW.m-2 
gDW.m'^ 
g  DW.m‘^ 
gDW.m'^ 
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TWRTG 

TWSD1,2,3 

TWSG1,2,3 

TWSTD 

TWSTG 

WTMP 


Total  dry  weight  of  live  roots  2  or  3  cohorts 
Total  dry  weight  of  dead  stems  2  or  3  cohorts 
Total  live  stem  dry  weight  cohort  1,2,3 
Total  dry  weight  of  dead  stems  2  or  3  cohorts 
Total  dry  weight  of  live  stems  2  or  3  cohorts 
Daily  water  temperature 


g  DW.m'^ 
g  DW.m'^ 
g  DW.m'^ 
g  DW.m'^ 
g  DW.m'^ 
°C 


B2 


Appendix  B  Outpn*  Parameters  Available 
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Public  reporting  burden  for  this  collection  of  infontnation  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  the  collection  of  infonnation.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions 
for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Infonnation  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Adington,  VA  22202-4302,  and  to  the 
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